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■ We study the non-linear dynamics of an inflationary phase transition in 

a quartically self coupled inflaton model within the framework of a de Sitter 

JL background. Large N and Hartree non-perturbative approximations com- 

Uh, 

bined with non-equilibrium field theory methods are used to study the self- 
consistent time evolution including backreaction effects. We find that when 
the system cools down from an initial temperature T, > T c to below T c with 
the initial value of the zero mode of the inflaton <^>(0) << mX" 1 ^ 4 , the dy- 
namics is determined by the growth of long- wavelength quantum fluctuations. 
For 4>(0) >> toA _1//4 the dynamics is determined by the evolution of the clas- 
sical zero mode. In the regime where spinodal quantum fluctuations give the 
most important contribution to the non-equilibrium dynamics, we find that 
they modify the equation of state providing a graceful exit from the infla- 
tionary stage. Inflation ends through this new mechanism at a time scale 
t s ^ [H/m 2 ] In [A -1 ] which for H > m and very weak coupling allows over 
one hundred e-folds during the de Sitter phase. Spatially correlated domains 
grow to be of horizon-size and quantum fluctuations "freeze-out" for times 
t > t s . 
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I. INTRODUCTION 



For all of its observational success, the standard big bang cosmology is not without 
its shortcomings. The inflationary scenario was conceived |I| to overcome many of these 
shortcomings, chief amongst them the flatness and horizon problems. The original old 
inflationary scenario |l| required a first-order phase transition with strong supercooling, 
during which the scale factor grew exponentially thus solving the homogeneity and flatness 
problem. This model, however, lacked any mechanism to provide a succesful exit from the 
inflationary phase B. The new inflationary or 'slow-rollover' scenario |],[|] overcame these 
problems. In these scenarios, a phase transition from an unbroken phase to a symmetry 
broken one takes place, and the scalar field, whose expectation value serves as the order 
parameter, drives inflation. This scalar field, the inflaton, drives a stage of exponential (de 
Sitter) expansion while it slowly rolls down the potential hill. Other versions, such as chaotic 
inflation || in which a phase transition is not necessary have been advocated (for reviews 
of inflationary cosmology, see @-[|). 

Phase transitions play an important role in the early universe, and, as mentioned above, 
are pivotal to many inflationary models. They are also necessary in models in which topo- 
logical defects provide the seeds for density perturbations |]10[ . It has been pointed out 
by many authors that using flat space-time effective potentials and equilibrium concepts to 
describe the dynamics of the inflationary phase transition can be misleading. Both gravita- 
tional effects as well as those due to quantum and thermal fluctuations have been studied 
and found to be very important for the description of the dynamics flTT|-|T6 . 



Although there have been previous studies of the dynamics of phase transitions JTl|-|To|, 
most of the attempts were hampered by the lack of a self-consistent non-perturbative treat- 
ment of the non-linearities and backreaction effects. It is only recently that non-equilibrium 
methods have been developed to study real time phenomena during phase transitions [|T7|-|2~0| 
and other strongly out of equilibrium situations pH] , including non-perturbative treatments 
22|^4|. During the last few years these methods have been adapted to study strongly 



out of equilibrium phenomena in cosmology, allowing the understanding of the dynamics of 



non-linear field theory in such situations p5| , |26 



Studies of the dynamics of typical second order phase transitions in Minkowski space-time 
have revealed a wealth of new and interesting phenomena. In situations in which the phase 
transition occurs on time scales which are faster than the relaxation time of long- wavelength 
fluctuations, the dynamics is driven by the growth of long-wavelength, spinodal instabili- 
ties. The backreaction of these quantum and thermal fluctuations dramatically modifies 
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the evolution of the zero mode of the scalar field eventually shutting off the instabilities 
|27j p9|. These instabilities are the hallmark of the phase transition and are responsible for 
the formation and growth of correlated domains [ZR2j|. 

The nature of these instabilities and their intrinsic non-perturbative dynamics is rather 
simple to understand p7|-|29|j . As the system cools down from a quasi-equilibrium disordered 
high temperature phase to below the critical temperature the effective mass of the scalar 
field becomes negative. If the expectation value of the field was zero in the equilibrium 
high temperature phase, it will remain zero below the critical temperature in the absence 
of biasing fields. The fluctuations of the field will grow in time and eventually become large 
enough to sample the minima of the (effective) potential, which are non-perturbatively large 
in amplitude. This growth of fluctuations translates into a growth of spatial correlations 
of the field, resulting in domains p8 |. In Minkowski space-time there is no restriction to 
the final size of these domains, and generically in weakly coupled theories these correlated 



regions reach sizes several times larger than the zero temperature correlation length [28 



In a de Sitter spacetime or in a general FRW cosmology, for that matter, causality limits 
the physical distance within which the order parameter can be correlated; in a de Sitter 
spacetime, this distance is dn = H~ x , with H the Hubble constant. 

The goal of this article is to study the full non-equilibrium dynamics of an inflationary 
phase transition including the backreaction effects in a simple model of a scalar field (the 
inflaton) with a quartic self-interaction. We use the methods of non-equilibrium field theory 
[T7| combined with self-consistent, non-perturbative Hartree and large iV approximations 

2i],miii,niis 

The metric is taken to be a fixed cosmological background within which we study the 
effects of quantum and thermal fluctuations on the evolution of the expectation value of the 
inflaton and the energy momentum tensor, both analytically and numerically, so as to follow 
the equation of state. When the temperature falls below the critical value, long- wavelength 
fluctuations grow exponentially as a result of the spinodal instabilities and their contribution 
to the equations of motion eventually become of the same order as the tree-level terms. We 
argue that in the very weakly coupled case, A ~ 10~ 12 , the change in temperature and the 
effective mass occur on time scales much shorter than the time in which long-wavelength 
fluctuations can adjust to local thermodynamic equilibrium, and they fall out of equilibrium 
just as in a "quench" from the high temperature phase. 

When these fluctuations are incorporated self-consistently as a backreaction in the evo- 
lution equations, we find, for slow-roll initial conditions, that the effects of these fluctuations 
dramatically change the dynamics. The unstable quantum fluctuations modify the equation 
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of state away from the vacuum dominated state required for de Sitter expansion, and for a 
reasonably wide range of parameters and initial conditions, these fluctuations provide a new 
mechanism to gracefully exit the inflationary stage within an acceptable number of e-folds. 
During this time, spatial correlations grow to reach horizon-size and freeze-out at larger 
times. 

In the large N approximation, we find that the late time dynamics is summarized by a 
very useful sum rule relating the infinite time limit of the order parameter 0(oo) and the 
infinite time limit of the quantum fluctuations (ip 2 (oo)) as follows, 

-|m 2 | + 2 (oo) + (^ 2 (oo)) = . 

This result implies that the excitations are massless and minimally coupled. However we 
find the new result that despite the fact that the relevant degrees of freedom are massless 
and minimally coupled asymptotically, their equal time two-point function saturates to a 
constant value rather than growing linearly in time. 

In section II we introduce the model, the approximations, and analyze the subtle but 
important issues of renormalization. Section III is devoted to presenting the renormalized 
equations of motion and the correlation functions. Section IV and V provide an analytic and 
numerical study of different well as an estimate of the backreaction of the quantum 

fluctuations onto the scale factor. Section VI presents a discussion of the results and the 
final conclusions and further avenues of study. 



II. THE MODEL AND THE APPROXIMATIONS 



The need for a non-perturbative treatment of nonequilibrium quantum field dynamics 
stems from the fact that when the temperature falls below the critical point, the effective 
mass term becomes negative. This results in an instability of the long-wavelength modes 
which grow exponentially for early times after the transition. The fluctuations of the scalar 
order parameter will grow in time to sample the broken symmetry vacua, leading to large 
amplitude fluctuations that can only be treated within a non-perturbative scheme. 

There are two approximation schemes that have been used to study the non-equilibrium 
dynamics during phase transitions, each with its own advantages and disadvantages. The 
Hartree factorization [|rj],^3],|2B],^S| has the advantage that it can treat the dynamics of a 
scalar order parameter with discrete symmetry, while its disadvantage is that it is difficult 
to implement consistently beyond the lowest (mean field) level. The advantage of the large N 
approximation pl],|22|,|25|,|26|,|30| is that it allows a consistent expansion in a small parameter 
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(1/iV) and correctly treats continuous symmetries in the sense that it implements Goldstone's 
theorem. Moreover, the large N expansion becomes a X/N expansion for small values of A. 
Therefore, it may be a reliable approximation for the typical values of A used for inflation 
even when N — 1. It should be noted that when spontaneous symmetry breaking is present, 
the large N limit always produces massless Goldstone bosons. 

Both methods implement a resummation of a select set of diagrams to all orders and 
lead to a system of equations that is energy conserving in Minkowski space time, and as will 
be shown below, satisfies covariant conservation of the energy momentum tensor in FRW 
cosmologies. Furthermore, both methods are renormalizable and numerically implementable. 
Given that both methods have advantages and disadvantages and that choosing a particular 
scheme will undoubtedly lead to criticism and questions about their reliability, we will study 
both, comparing the results to obtain some universal features of the dynamics. 

We restrict our study to a spatially flat Friedmann- Robertson- Walker universe with scale 
factor, a(t), and metric, 

ds 2 = dt 2 - a 2 (t) dx 2 , (2.1) 
and focus in particular on the case of a de Sitter space-time with scale factor, 

a(t) = e m . (2.2) 
The action and Lagrangian density are given by, 

S = Jd 4 x C, (2.3) 



C = a 3 (t) 



1=>\ . l(V$(a;)) 2 

-$ (x) V($(x)) 

2 V ' 2 a 2 (t) y y n 



(2.4) 



t ,^« A /-j 2 2NM 2 
K J 8N \ A 



M 2 = -m 2 +£K, 



(2.5) 



where we have included the coupling of to the scalar curvature TZ(t) since it will arise as 
a consequence of renormalization. In the de Sitter universe ( |2.2j ), 1Z = 12H 2 . The canonical 
momentum conjugate to $(x) is, 
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ri( x ) = a 3 (t) (2.7) 
and the time dependent Hamiltonian is given by, 

H{t) = I dix {^) + ^^ x))2 + a3{t)V ^}- (2 ' 8) 



A. The Hartree Approximation 

To implement the Hartree approximation, we set N = 1 and write, 

${x,t) =<f){t)+ij;{x,t), (2.9) 

with, 

0(t) = <$(£,*)>; =0, (2.10) 

where the expectation value is defined by the non-equilibrium density matrix specified below, 
and we have assumed spatial translational invariance, compatible with a spatially flat metric. 
The Hartree approximation is obtained after the factorization, 

V 3 (x,t)-,3(V 2 (x,t))V(x,t), (2.11) 

^ 4 (f,t)^6^ 2 (f,t))^ 2 (f,t)-3(^ 2 (f,t)) 2 , (2.12) 

where by translational invariance, the expectation values only depend on time. In this 
approximation, the Hamiltonian becomes quadratic at the expense of a self-consistent con- 
dition. 

At this stage we must specify the non-equilibrium state in which we compute the expec- 
tation values above. In non-equilibrium field theory, the important ingredient is the time 
evolution of the density matrix p(t) (see |T^] and references therein). This density matrix 
obeys the quantum Liouville equation whose solution only requires an initial condition p(tj) 



|l7| , pl|j20| , |2qj26| . The choice of initial conditions for this density matrix is an issue that 
pervades any calculation in cosmology. Since we want to study the dynamics of the phase 
transition, it is natural to consider initial conditions that describe the instantaneous modes 
of the time dependent Hamiltonian as being initially in local thermodynamic equilibrium 
at some temperature Tj > T c . Given this initial density matrix, we then evolve it in time 
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using the time dependent Hamiltonian as in p5] or alternatively using the complex time 



path integral method as described in [ p^ , p^ - pl| , p6 

Following the steps of references pj|^? 
tion value of the inflaton field to be, 

A 



we find the equation of motion for the expecta- 



3A 



<f>(t) + 3Hc/>(t) + M 2 </>(t) + ^\t) + — 0(t){^ 2 (t)> 



0. 



(2.13) 



The equal time correlation function is obtained from the coincidence limit of the non- 
equilibrium Green's functions, which are obtained from the mode functions obeying 



fk(t) = 0, 



(2.14) 



with the effective frequencies, 



a 2 (t) 



+ M 2 (t) , 



where 



M 2 (t) 



M 2 + y0 2 (t) + y<^(t)>. 



(2.15) 



(2.16) 



At this stage we must provide the initial conditions on the mode functions fk(t). As 
mentioned above our choice of initial conditions on the density matrix is that of local ther- 
modynamic equilibrium for the instantaneous modes of the time dependent Hamiltonian at 
the initial time. Therefore we choose the initial conditions on the mode functions to rep- 
resent positive energy particle states of the instantaneous Hamiltonian at t — 0, which we 
chose as the initial time. Therefore our choice of boundary conditions at t = is 

1 



A(o) 



/*(o) 



W k ; Wk 



k 2 + Ml 



o • 



(2.17) 



where the mass Mq determines the frequencies c^fc(0) and will be obtained explicitly later. 
With these boundary conditions, the mode functions /fc(0) correspond to positive frequency 
modes (particles) of the instantaneous quadratic Hamiltonian for oscillators of mass M . 
The equal time correlation function is given in terms of the mode functions as 
91, 



21,25 



d 3 k \f h (t)f 



(2tt) 3 2 

The energy and pressure density are given by pB] 



coth 



Wk 

2T, ; 



(2.18) 
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1 -,2/ x A / ,0/ \ 2M 



2\ 2 



+ 



2(2tt) : 



coth 



22} 



\f k (t)\ 2 + u 2 k (t)\f k (t)f 



3A 



(2.19) 



p + e = 2 (t) + 



2(2vr) 



coth 



22} 



|A(0I 2 + 



3a 2 (t) 



(2.20) 



It is straightforward to show using the equations of motion ( |2.13| , p7T^ ) that the bare 
energy is covariantly conserved, 



e + 3H(p + e) = 0. 



(2.21) 



B. The Large N Approximation 

To obtain the proper large N limit, the vector field is written as, 

$(£,£) = (a(x,t),7f(x,t)), 

with 7? an - 1-plet, and we write, 

a(x,t) = VN<f>(t)+ X {x,t) ; (a(x, t)) = VN<f>{t) ; ( x {x,t)) = 0. (2.22) 

To implement the large N limit in a consistent manner, one may introduce an auxiliary 
field as in However, the leading order contribution can be obtained equivalently by 

invoking the factorization, 



x 4 - 


-» 6(x 2 )x 2 + constant, 


(2.23) 






(2.24) 


7T • 7TJ - 


^2(7r 2 )7r 2 - (7r 2 ) 2 + 0(l/AT), 


(2.25) 


->2 2 
7T X ~ 




(2.26) 


T? 2 X ~ 


- (v? 2 )x- 


(2.27) 



To obtain a large N limit, we define, 

N-l 

7f(x,t)=i/j(x,t) (1,1, •••,!), (2.28) 
where the large A/" limit is implemented by the requirement that, 
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(2.29) 



The leading contribution is obtained by neglecting the 0(1/ N) terms in the formal limit. 
The resulting Lagrangian density is quadratic, with a linear term in x- The equations of 
motion in this case become, 



4(t) + 3H<P(t) + M 2 0(t) + ^(t) + \<P{t)(4> 2 {t)) 



0. 



(2.30) 



with the mode functions, 



d 3 k \f k {t)\< 
(2tt) 3 2 



coth 



Wk 

m 



3H 



dt 



<4V) 



fk(t) = o, 



and the effective frequencies, 



a 2 (t) 



+ M 2 {t) , 



where, 



M\t) = M 2 + ^(t) + ^(t)). 



(2.31) 



(2.32) 



(2.33) 



(2.34) 



The initial conditions are chosen to reflect the same physical situation as in the Hartree 
case, that is, the instantaneous particle states of the Hamiltonian at t — are in local 
thermodynamic equilibrium at some initial temperature higher than the critical value. Thus, 
as in the Hartree case but with modified frequencies, the initial conditions at t — are chosen 
to describe the instantaneous positive energy states, 

1 



A(o) 



; A(0) = -i ^W k ; W k = y/k* + Ml 



[2.35) 



We have maintained the same names for the mode functions and Mq to avoid cluttering 
of notation; their meaning for each case should be clear from the context. Notice that 
the difference between the Hartree and large iV case is rather minor. The most significant 
difference is that, in the equations for the zero modes, the Hartree case displays a factor 3 
difference between the tree level non-linear term and the contribution from the fluctuation 
as compared to the corresponding terms in the large iV case. The equations for the mode 
functions are the same upon a trivial rescaling of the coupling constant by a factor 3. 

The particular case with (f>(t) = is of interest, since by symmetry, it is a fixed point of 
the dynamics of the zero mode and corresponds to the case in which the phase transition 
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occurs from the symmetric phase into the broken phase in the absence of symmetry breaking 
perturbations or initial bias in the field. In this important case, the Hartree and large N 
mode equations are identical after the coupling constant is rescaled by a factor of 3, and 
our conclusions will be universal in the sense that both the large iV approximation and the 
Hartree approximation describe the same non-perturbative dynamics. In the large N limit 
we find the energy density to be given by, 



6 1 12^ A / j2 /x 2M ^ 



+ 



d 3 k 

2(2vr) : 



coth 



Wk 

2T,; 



l/*(f)l 2 +"2(f)l/*(*)l 5 



A 



(2.36) 



and (p + e)/N has the same form as for the Hartree case above (eq. ( |2.20| )), but in terms 
of the mode functions obeying the large N equations. Again, it is straightforward to show 
that the bare energy is covariantly conserved by using the equations of motion for the zero 
mode and the mode functions. Note that equation ([2.21 ) also holds in the large N limit. 

A relation analogous to ( |2.21| ) holds for the integrands of e and e + p in eqs. ( [2.1 9| ) and 
( CT ). Namely, 

d 



s k (t) + 3Ha k (t) = -[M 2 (t)] \f k (t)\ 



[2.37) 



Here. 



e k (t) = \f k (t)\ 2 + u;l(t)\f k (t)\ 

k 2 



a k (t) = 2 



!/*(*) r + 



3a 2 (t) 



\fkW 



(2.38) 



Equation ( |2.37| ) will be used below to prove that the renormalized energy is conserved. 



III. RENORMALIZED EQUATIONS OF MOTION 



Renormalization is a very subtle but important issue in gravitational backgrounds 
13| , |32|j33| . The fluctuation contribution (ip 2 (x,t)) , the energy and the pressure all need 



to be renormalized. The renormalization aspects in curved space times have been discussed 
at length in the literature |T3|,|3"2"|,|3"3"| and has been extended to the Hartree and large N 



self-consistent approximations for the non-equilibrium backreaction problem in f25| , |30 



In terms of the effective mass term for the large N limit given by ( |2.34|) and the Hartree 
case, eq. ( |2.16| ), and defining the quantity, 
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B(t) = a 2 (t) [M 2 (t)-TZ/6 



(3.1) 



we find the following large k behavior for the case of an arbitrary scale factor a(t) (with 
a(0) = 1), 

1 1 



\fkW 



+ 



ka 2 (t) + 2k 3 a 2 (t) 
1 



H 2 (0) - B{t) 



d 



|A(*)I S 



la(t) 2 k 5 
k 



B(t)[3B(t) - 2H 2 (0)} + a(t)- [a(t)B(t)\ + D + £>(l/£; 7 ) 



a 4 (t) 



+ 



1 



+ 



fca 2 (t) 
{-B(f) 



a(tyB(t) + 3a(t)a(t)B(t) - 4a 2 (t)B(t) 



3a(t) 4 A; 3 

+ 2# 2 (0)[2d 2 (t) + £(t)] + D } + G(l/A; 5 ). (3.2) 

The constant Do depends on the initial conditions and is unimportant for our analysis. 

Although the divergences can be dealt with by dimensional regularization, this procedure 
is not well suited to numerical analysis. We will make our subtractions using an ultraviolet 
cutoff constant in physical coordinates. This guarantees that the counterterms will be time 
independent. The renormalization then proceeds much in the same manner as in reference 
E3|; the quadratic divergences renormalize the mass and the logarithmic terms renormalize 



coupling constant and the coupling to the Ricci scalar. The logarithmic subtractions can be 
neglected because of the coupling A ~ 10~ 12 . Using the Planck scale as the cutoff and the 
inflaton mass tur as a renormalization point, these terms are of order Xhi[M p i/mR] < 10~ 10 , 
for m > 10 9 GeV. An equivalent statement is that for these values of the coupling and 
inflaton masses, the Landau pole is beyond the physical cutoff M v \. Our relative error in 
the numerical analysis is of order 10 -8 , therefore our numerical study is insensitive to the 
logarithmic corrections. Though these corrections are fundamentally important, numerically 
they can be neglected. Therefore, in what follows, we will neglect logarithmic renormaliza- 
tion and subtract only quartic and quadratic divergences in the energy and pressure, and 
quadratic divergences in the fluctuation contribution. 

Using the large k behavior of the mode functions, we find the finite temperature factor 
to be 127 



d 3 k- 



1 

24 



T? 



a 2 (t) 



[1 + O(M /T t ) + O (MMo/^)) + ■ 



(3.3) 



where we have used Tj > T c >> M . The first term in ( |3.3j ) arises from the l/[ka 2 (t)} in the 
asymptotic limit (p.2|). For the large iV case, we find the renormalization condition, 
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A 



- < + '-f<f(t) + £bK + '-fmt)) B = mUt, Ti) + '-f<p 2 (t) + t R 7L + '-f(^{t)) Ri (3.4) 
with the effective time dependent mass term in de Sitter space time, and subtracted fluctu- 



A 



A 



R 



ation contribution [25 



m 2 R (t,Ti) 



m 



R 



J^2 



e -2Ht _ 1 



R 



4tt 2 



k 2 dk coth 



2T, 



{Mt)\ l 



ka 2 (t) 2k 3 a 2 (t) 



H 2 (0)-a 2 (t) (M 2 (t)-TZ/Q 



(3.5) 



where k is a renormalization point, which for convenience will be chosen as k = \m R \, i.e. 
the renormalized inflaton mass. The Hartree case is similar but with a rescaling of the 
coupling constant by a factor 3. We will fix the renormalized coupling to the Ricci scalar to 
be = 0, thus ensuring minimal coupling in the renormalized theory. 

Assuming that Tj/T c w 0(1), we see that the phase transition occurs within the first 
few e-folds of inflation. When the temperature falls below the critical value, the effective 
mass becomes negative. As will be seen explicitly below, when this occurs, long- wavelength 
modes become unstable and grow. Local thermodynamic equilibrium will set in again if 
the contribution from the quantum fluctuations can grow and adjust to compensate for 
the negative mass terms on the same time scales as that in which the temperature drops. 
However, as discussed below, for very weak coupling the important time scales for the non- 
equilibrium fluctuations are of the order of [H/m 2 ^] ln[l/A], which are much longer than the 
time it takes for the temperature to drop well below the critical value to practically zero. 
Thus, the non-equilibrium dynamics will proceed as if the phase transition occured via a 
"quench" , that is with an effective mass term, 



rn 



eff 



(t) = m 2 e{t L -t)- m\Q{t - 1,) 



m 



R 



c 



- 1 



Therefore we choose the initial conditions on the mode functions at 
terms of the effective mass, 



Ml 



m 



R 



T 2 

- 1 

^2 



A 



^0 2 (O) ; §f >i, 



>0. (3.6) 
to be given in 

(3.7) 



with X R —>■ 3Ar for the Hartree case. 

We introduce the following dimensionless quantities and definitions, 

H k 

m R ' q m R ' 



T 



m R t 



h 



(3.8) 



12 



Ti Ti W k A 

r = — ; T = — ; u q = — ; g 



T, 



m R 



m R 



8tt 2 ' 



/,(T) = v^A(t) I A' = 



ft 



With £_r = and X = 1 the equations of motion become: 
Large N: 

ri + 3hf] — r] + r] 3 + (?E(r)?7 = 



d 2 3h d 
dr 2 dr 



a 2 (r) 



l + V 2 + gE(r) 



f q (r) = 



(3.9) 



(3.10) 



(3.11) 



(3.12) 
(3.13) 



Hartree: 



fOJn 



; /,(0) = -i v ^ ; o;, = ^2 + 7 .2_ 1+r7 2(o) 



57 + 3/1?) — 77 + r? 3 + 3(?S(r)r] = 



d 2 o , d g 2 
+ 3/i— + 



/ 9 (0) 



dr 2 dr o 2 (t) 
1 



1 + 3r] 2 + 3#E(r) 



/ ff (r) = 



/UJn 



; f q $) = -iy/UT q ; ^ = ^ + ^-1 + 3^(0) 



The dots stand for derivatives with respect to r, and in both cases, 



0E(r) 



g 2 dg coth 



2T 



2g 3 a 2 (t) 



^ 2 (0) 



a 2 {t) 



m 



R 



(M 2 (t)-K/Q) 



(3.14) 

(3.15) 
(3.16) 

(3.17) 



. (3.18) 



Numerically, the most significant contribution to (ip 2 ) arises from low wavevectors q < 
10 — 20 in all of the cases studied (see figures l.e and l.f below). Since > T c = ^2Am 2 R /\, 
q « T and M << Tj for these low momentum modes, they are "classical" and we can 
approximate, 



coth 





2T 


2r I 


12T. 




u q \ 



(3.19) 
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In this approximation, the fluctuation contribution becomes, 



9%(t) = g 



q 2 dq 



124 



\f q (r)\" 



2q 3 a 2 (t) 



a 2 (t) 



h 2 (Q)- ^(M 2 (t)-n/Q 



m 



R 



g = 2rg\ — = 1.103 r^fg = 0.124 rVX . 
V a 

The renormalized energy density and pressure for the large N case are given by, 



(3.20) 



e 2m|j j V 2 1 



N A [24 

% - i) m 

q 3 



+ 7(^-ir+ 



2 , 9 fQ 



\Ur)\ 2 + ^ 2 Ar)\f q {r)f 



2q a(t) 



a(t) 4 



E(r) 



(3.21) 



p + e 
N 



with, 



2m% 



f] +g 



q 2 dq 



\f q (r)\ 2 + 



q2 \m\ 2 - Aq 



l{t) 0(q - 1) 



3a 2 ( 



T 



3a(ty 



.2/ N 1 , , „2 



<(r) 



a 2 (r) 



-l + ^(r)+^S(r) 



5(0 

(3.22) 
(3.23) 



and B{t) is given by eq. ( |3.1| ). The coefficients a(t), f3(t), 7(i) and <5(t) are obtained from 
the asymptotic behaviour of Skit) and 0>(£) (|2.38| - |3~2|) . We find, 

«(*) = 2 1 /,U + + ^(°) 2 > 

m R a{t) 

M) = 1 4 1 m4 + 2a(f)a(f)fl(f) 

4m|ja(t) 4 L 

- 2[a 2 (t) + a 2 (£) M(£) 2 ][B(i) - #(0) 2 ] + A)} , 

-l2^F{ a2(t,fi(i, - 5a(t) " (i,B(i) 

+ 6d 2 (t) [B(t) - H(0) 2 } -2H(0) 2 B(t) - 2D ] . 



5(t) 



(3.24) 



For the Hartree case (setting N = 1) the only changes are that we take g£/4 — > 3gE/4 
in the last term in the energy, and use the frequencies, 
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^ 2 ( r ) = ~TTT " 1 + 3v 2 (t) + 3gX(r), (3.25) 



and that the initial conditions on the mode functions are given by (|3.17| ). 

In order to prove that the renormalized energy and pressure fullfils the continuity equa- 
tion ( |2.21| ), we need to study the properties of the subtracted terms in equations ( |3"3 



( |3.22| ). Inserting the asymptotic behaviour of Ek(t), auit) and \fk(j)\ 2 for large k in ( p.37|) , 
we find, 

a( t ) + 6// 7W = ^|[M=( t )], 
3(t) + 6H S(t) = ^ \H(0f - B(t)]j t [Af\t)} . (3.26) 



Using these relations and (|2.37 ), it is straightforward to show that the renormalized energy 



and pressure given by eqs. (|3.21|) - (|3.22|) satisfy the continuity equation, 

e + 3H(p + e) = . 

A noteworthy point is that when the cutoff is kept fixed in physical coordinates, upon 
taking the time derivative in the integrals, there is a contribution from the upper limit of 
the integrals. However the subtractions guarantee that in the formal limit when the cutoff is 
taken to infinity this contribution vanishes. While the existence of the Landau pole beyond 
the Planck mass restricts taking the cutoff to infinity in this effective theory, we find that 
any contributions from the upper limit are numerically small. In our numerical evolution, 
the energy density is covariantly conserved to one part in 10 7 . 

From the evolution of the mode functions that determine the quantum fluctuations, we 
can study the growth of correlated domains with the equal time correlation function, 

S(x,t) = «.(f,*)«(0,i)>, (3.27) 

(2tt) 3 2 \2Ti, 

which can be written in terms of the power spectrum of quantum fluctuations, 

5 ( M ) = £fal) ; 5( ff ,r) = |/„(r)| a ooth(^:). (3.28) 
m R V2T / 

It is convenient to define the dimensionless correlation function, 

S(\x\,t) 1 r 

S(p,r) = ^— = —- / qdqsm[qp]S(q,r) ;p = m R \x\. (3.29) 

m R 4tt z p J 

We now have all the ingredients to study the particular cases of interest. 
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IV. EVOLUTION FOR 0(0) = 0(0) = 



A. Analytical Results 



We begin by considering the situation in which the expectation value of the inflaton field 
sits atop the potential hill with zero initial velocity. This situation is expected to arise if 
the system is initially in local thermodynamic equilibrium an initial temperature larger than 
the critical temperature and cools down through the critical temperature in the absence of 
an external field or bias. 

The order parameter and its time derivative vanish in the local equilibrium high tem- 
perature phase, and this condition is a fixed point of the evolution equation for the zero 
mode of the inflaton. There is no rolling of the inflaton zero mode in this case, although the 
fluctuations will grow and will be responsible for the dynamics. 

We can understand the early stages of the dynamics analytically as follows. For very 
weak coupling and early time we can neglect the backreation in the mode equations, which 
in both the large N and Hartree cases become, 

^ + 3 ^ + ^)- 1 ]^-°' ^ 



/,(0) 



; f q (r) 



The solutions are of the form, 



q 2 + r 2 — 1 . 



1 



3 q 
f q (r) = exp[--hr] {a(q) J v (z) + b(q) J-„(z)} ; z = -exp[-/ir] ; v= 2 

where the coefficients a(q) and b(q) are determined by the initial conditions: 



b(q) 



irq 



2h sin z/7T 



iu q — \ h 



9 
4' 



(4.2) 



(4.3) 



(4.4) 



a(q) 



7r q 



2h sin z/7T 



™*--* h J-A%\-JL 



1 

v \h 



q \hy 

For long times, e hr > q/h, these mode functions grow exponentially, 

6(g) (2h 



f q (r) ~ b(q) J_ v (z) 



(4.5) 



(4.6) 



V{l-v) V q J 

The Bessel functions appearing in the expression for the modes f q (j) can be approxi 
mated by their series expansion, 
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This is an expansion in powers of q 2 /{yh 2 ) and we conclude that gSfV) is dominated by the 
modes with q < \fh. 

The integral for (?S(r) can be approximated by keeping only the modes q < fVh, 
where / is a number of order one, and by neglecting the subtraction term which will cancel 
the contributions from high momenta. Numerically, even with the backreaction taken into 
account, the integral is dominated by modes q < f ~ 10 — 20 in all of the cases that we 
studied (see fig. (l.e) below). 

The contribution to the fluctuations from these unstable modes is: 

9 V ' ~ V 6 4tt 2 Mi \ ml ) ' y ' 

where again, we have taken the high temperature limit, Tj ~ T c ^> itlr. 

From this equation, we can estimate the value of r s , the "spinodal time", at which the 
contribution of the quantum fluctuations becomes comparable to the contribution from the 
tree level terms in the equations of motion. This time scale is obtained from the condition 
jS(t.) = 0(1): 



1 In 



(2u-3)h 



f 3 h 3/2 T 1 ( + M 2 



6 4tt 2 m R Ml T c \ m% 



(4.9) 



which is in good agreement with our numerical results, as will become clear below (see 
figures (l.a), (2. a) and (3.b)). For values of h > 1, which, as argued below, lead to the most 
interesting case, an estimate for the spinodal time is, 

r 8 ~yln[l/VA] + 0(l) (4.10) 

which is consistent with our numerical results (see figure (l.a)). 

For r > t s , the effects of backreaction become very important, and the contribution from 
the quantum fluctuations competes with the tree level terms in the equations of motion, 
shutting-off the instabilities. Beyond r s , only a full numerical analysis will capture the 
correct dynamics. 

It is worth mentioning that had we chosen zero temperature initial conditions, then the 
coupling g — > g (see (|3.20|) ) and the estimate for the spinodal time would have been, 

r s ~^m[l/A] + C>(l), (4.11) 



17 



that is, roughly a factor 2 larger than the estimate for which the de Sitter stage began at 
a temperature above the critical value. Therefore ( |4.10| ) represents an underestimate of the 
spinodal time scale at which fluctuations become comparable to tree level contributions. 

The number of e-folds ocurring during the stage of growth of spinodal fluctuations is 
therefore, 



QU2 

JVe»— ln[l/VA], (4.12) 



or in the zero temperature case, 



3/? 2 

ln[l/A], (4.13) 



which is a factor 2 larger. Thus, it becomes clear that with A ~ 10~ 12 and h > 2, a required 
number of e-folds, M e ~ 100 can easily be accomodated before the fluctuations become large, 
modifying the dynamics and the equation of state. 

The implications of these estimates are important. The first conclusion drawn from these 
estimates is that a "quench" approximation is well justified (see figure (l.a)). While the tem- 
perature drops from an initial value of a few times the critical temperature to below critical 
in just a few e-folds, the contribution of the quantum fluctuations needs a large number of 
e-folds to grow to compensate for the tree-level terms and overcome the instabilities. Only 
for a strongly coupled theory is the time scale for the quantum fluctuations to grow short 
enough to restore local thermodynamic equilibrium during the transition. 

The second conclusion is that most of the growth of spinodal fluctuations occurs during 
the inflationary stage, and with A ~ 10~ 12 and H > m#, the quantum fluctuations become 
of the order of the tree-level contributions to the equations of motion within the number of 
e-folds necessary to solve the horizon and flatness problems. Since the fluctuations grow to 
become of the order of the tree level contributions at times of the order of this time scale, 
for larger times they will modify the equation of state substantially and will be shown to 
provide a graceful exit from the inflationary phase within an acceptable number of e-folds. 

For r < t s , when the contribution from the renormalized quantum fluctuations can be 
ignored, the Hubble constant is given by the classical contribution to the energy density. In 



terms of the dimensionless quantities introduced above ( |3.10| ), we have, 

167rml 



_ OR 

3AM 2 |_ 



(4-14) 



In the situration we consider here, with rj = rj = 0, the condition that h > 2 for A ~ 10 12 
translates into ~ 10 13 GeV , which is an acceptable bound on the inflaton mass. 
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To understand more clearly whether or not the effect of quantum fluctuations and growth 
of unstable modes during the inflationary phase transition can provide a graceful exit sce- 
nario, we must study in detail the contribution to the energy and the equation of state of 
these quantum fluctuations. 

Although we are working in a fixed de Sitter background, the energy and pressure will 
evolve dynamically. A measure of the backreaction effects of quantum fluctuations on the 
dynamics of the scale factor is obtained from defining the 'effective Hubble constant', 



Therefore, the quantities, 



n(r) 

W(0) \e(0)' 



(4.16) 



and 



H{t) _3 
H\r) ~ 2 



1 + P(T) 



(4.17) 



s(r)_ 

give dynamical information of the effects of the backreaction of the quantum fluctuations 
on the dynamics of the scale factor. Whenever p(r) + e{r) ^ 0, 7Y(r)/7Y(0) 7^ 1, or 
H(r)/7i 2 (r) 7^ 0, the backreaction from the quantum fluctuations will dramatically change 
the dynamics of the scale factor, and it will no longer be consistent to treat the scale factor 
as fixed. When 7i(r) /7i(0) <C 1 the de Sitter era will end. 

From this point onwards only a full treatement of the backreaction, including the correct 
dynamics of the scale factor, will describe the physics. The time scale on which the quantum 
fluctuations will begin to influence the dynamics of the scale factor is of the order of the 
spinodal time estimated above, since the contribution of the quantum fluctuations becomes 
comparable to the tree level terms and modifies the equation of state. 

Therefore, there is the possibility that the growth of quantum fluctuations can provide 
a graceful exit from the inflationary phase, even when the zero mode does not roll. The 
parameters should be chosen in such a way so that the requisite 60 or more e-folds of 
expansion take place before the spinodal time. From the estimates provided above, this is 
relatively easy to accommodate with reasonable values of the inflaton mass and for the weak 
coupling that is usually assumed in inflationary models. 
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B. Numerical Analysis 



We now solve the large N set of equations ( |3.13|) numerically, with the initial conditions 
([TLJ), taking 77(0) = 77(0) = 0. 

The numerical code is based on a fourth order Runge-Kutta algorithm for the differential 
equation and an 11-points Newton-Cotes algorithm for the integral, with a typical relative 
errors 10~ 9 in the differential equation and in the integrals. We have tested for cutoff 
insensitivity with cutoffs q m ax ~ 50 ; 100 ; 150 with no appreciable variation in the numerical 
results. The reason for this cutoff insensitivity is due to the fact that only long- wavelength 
modes grow in amplitude to become non-perturbatively large, whereas the short-wavelength 
modes always have perturbatively small amplitudes. We have chosen r = Tj/T c = 2 as a 
representative value and A = 10~ 12 . The insensitivity on the value of the cutoff confirms 
that the high-temperature limit ( |3.19| ) is warranted. 

As argued previously, for A ~ 10~ 12 , the cosmologically interesting time scales for the 
spinodal instabilities to grow during say the first 60-100 e-folds of inflation occur for h > 1, 
leading to H > > 10 13 GeV which is a phenomenologically acceptable range for the 
Hubble constant during the inflationary stage. 

Figure (l.a) shows the contribution from the quantum fluctuations, (?£(t) vs. r for 
A = 10~ 12 ; Tj/T c = 2 ; h = 2 ; i](0) = ; ?)(0) = 0. The quantum fluctuations, as 
measured by gE(r), grow to be of order 1 in a time scale r w 40 which is the time scale 
predicted by the early time estimate ( f4.10|) . Figure (l.b) shows [p(r) + e(r)] A/(2mfj) vs. 
t for the same values of figure (l.a). Initially, p = —e and the quantity p + e is zero. 
At the spinodal time, there is a change in the equation of state, causing p + e to grow. 
However for late times, the energy density and pressure are each redshifted away such that 
the sum again approaches zero. We have checked numerically that the energy is covariantly 
conserved, obeying the relation e + 3H(p + e) = to our numerical accuracy of one part 
in 10 7 . Figures (l.c,d) show TC(r)/TC(0) and H.{t) /1-L 2 {t) vs. r respectively. These figures 
show clearly that when the spinodal quantum fluctuations become comparable to the tree 
level contribution to the equations of motion, the backreaction on the scale factor becomes 
fairly large. At this point, the approximation of keeping a fixed background breaks down 
and the full self-consistent dynamics will have to be studied. At this time, the inflationary 
stage basically ends since 7i is no longer constant. This occurs for r ~ 40 giving about 80 
e-folds of inflation during the time in which 7i is approximately constant and equal to H. 
Therefore, this new mechanism of spinodal fluctuations, with the zero mode sitting atop the 
potential hill provides a graceful exit of the inflationary era without any further assumptions 
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on the evolution of the scalar field. 

These fluctuations translate into an amplification of the power spectrum at long wave- 
lengths for q « h. To see this clearly we plot gS(q, r) with S(q, r) given by eq. ( ft.28| ) vs. 
q for t = 60 in fig. (l.e). This quantity is very small, because of the coupling constant in 
front, but for r « r s it grows to be of order one for long wavelengths (see also fig. (l.h)) and 
vanishes very fast for q > 10. The integral in gH{r) is dominated by these long wavelengths 
that become non-perturbatively large, whereas the contribution from the short wavelengths 
remains always perturbatively small. This is the justification for the approximations per- 
formed early that involved only the long-wavelength modes and cutoffs of order y/h. The 
equal time spatial correlation function given by eq. ( |3.29|) can now be computed explicitly. 



Figure (l.f) shows S(p; r)/S(0; r) as a function of p for r > 2. We define the correlation 
length £(t) as the value of p for which the ratio is 1/e. Figure (l.g) shows £(t); notice that 
the correlation length saturates to a value £(oo) wl/ft, and that the correlated regions are 
of horizon size. 

We have performed numerical analysis varying h with the same values of A and for the 
same initial conditions, and found that the only quantitative change is in the time scale 
for gTt(r) to be of order one. We find that the spinodal time scale grows almost linearly 
with h and its numerical value is accurately described by the estimate (|4.10|) . The case in 
which the Hubble constant is h = 0.1 is shown explicitly. Figure (2. a) shows g£(r), which 



demonstrates the oscillatory behavior similar to what is seen in Minkowski space |26[| . The 
correlation length, £(r), is shown in fig. (2.b); its asymptotic value is again approximately 
given by 1/h. 

C. The late time limit 

For times r > r s « 40 (for the values of the parameters used in figures (1)) we see from 
figures (l.a,h) that the dynamics freezes out. The fluctuation, ^S(r) = 1, and the mode 
functions effectively describe free, minimally coupled, massless particles. The sum rule, 

-l + #£(oo) = 0, (4.18) 



is obeyed exactly in the large N limit as in the Minkowski case [2B . 

For the Hartree case g — > 3g, but the physical phenomena are the same, with the only 
difference that the sum rule now becomes gE(oo) = 1/3. We now show that this value is 
a self-consistent solution of the equations of motion for the mode functions, and the only 
stationary solution for asymptotically long times. 
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In the late time limit, the effective time dependent mass term, — 1 + rf + gH, in the 
equation for the mode functions, ( |3.13| ), vanishes (in this case with rj = 0). Therefore, the 
mode equations ( |3.13| ) asymptotically become, 

J2 



d 2 | 3h d | r 

dr 2 dr a 2 (r) 



The general solutions are given by, 



asy. 



T 



cxp 



-hr 
2 



d(q) J3/2 ( je 



f q (r) = . 



+ c(q) N 3/2 I j-e 



(4.19) 



(4.20) 



where Jz/2{ z ) and N 3 / 2 { z ) are the Bessel and Neumann functions, respectively. The coef- 
ficients, d(q) and c(q), can be computed for large q by matching fg Sy (r) with the WKB 
approximation to the exact mode functions f q (r) that obey the initial conditions ([3.14 ). 
The WKB approximation to f q (r) has been computed in ref. f25(| , and we find for large q, 



d(q) 




-{h + A) + 0{q 
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e -iq/h + 
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e i( i/ h 
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(4.21) 



c(q) 
where 




-{h + A) + 0{q- 



e- iq/h + u 



Inh 
~8q 



l + O 



Jg/h 



poo 

A = / dr e hT M 2 (r) . 
Jo 

In the t — > 00 limit, we have for fixed q, 



3/2 



asy 1 



c(q) . 



(4.22) 



(4.23) 



(4.24) 



which are independent of time asymptotically, and explains why the power spectrum of 
quantum fluctuations freezes at times larger than the spinodal. This behavior is confirmed 
numerically: fig. (l.h) shows ln[|/ g (r)| 2 ] vs. r for q = 0,4,10. Clearly at early times the 
mode functions grow exponentially, and at times of the order of t 8 , when g^{r) ~ 1 the 
mode functions freeze-out and become independent of time. Notice that the largest q modes 
have grown the least, explaining why the integral is dominated by q < 10 — 20. 
For asymptotically large times, gll is given by, 



g-E(oo) = g h 2 f 
Jo 



* coth ( ^ 
q V2T 



2h 2 
— I c(q) 

7T 



q 



(4.25) 



where only one term in the UV substraction survived in the r = 00 limit. For consistency, 
this integral must converge and be equal to 1 as given by the sum rule. For this to be 
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the case and to avoid the potential infrared divergence in ( 4.25Q , the coefficients c(q) must 



vanish at q = 0. The mode functions are finite in the g — f limit provided, 

c{q) q =°Cq 3 / 2 , (4.26) 

where C is a constant. 

The numerical analysis and figure (l.e) clearly show that the mode functions remain 
finite as q — > 0, and the coefficient C can be read off from these figures. This is a remarkable 
result. It is well known that for free massless minimally coupled fields in de Sitter space- 
time with Bunch-Davies boundary conditions, the fluctuation contribution (if) 2 (x,t)) grows 
linearly in time as a consequence of the logarithmic divergence in the integrals |12|- 14 



However, in our case, although the asymptotic mode functions are free, the coefficients that 
multiply the Bessel functions of order 3/2 have all the information of the interaction and 
initial conditions and must lead to the consistency of the sum rule. Clearly the sum rule 
and the initial conditions for the mode functions prevent the coefficients d(q) and c(q) from 
describing the Bunch-Davies vacuum. These coefficients are completely determined by the 
initial conditions and the dynamics. This is the reason why the fluctuation freezes at long 
times unlike in the free case in which they grow linearly []T2"|-|I4|. 

It is easy to see from eqs. (|3.21|) - (|3.22|) and (|4.24j ) that the energy and presure vanish for 



r — > oo. 

Analogously, the two point correlation function can be computed in the late time regime 
using the asymptotic results obtained above. Inserting eq. (|4.20p for the mode functions in 
eq. (|3.29|) yields the asymptotic behavior: 



S(p,r) t =°° -1- f°°qdq singpcoth (£) ^ | c(q) | 2 . (4.27) 
Aix z p Jo \2TJ nq 6 

The asymptotic behavior in time of the equal time correlation function is thus solely a 
function of r. The large r behaviour of S(f, +oo) is determined by the singularities of 
| c(k) | 2 in the complex k plane. We find an exponential decrease, 

S(p,+oo) P ~ C , (4.28) 

P 

where p = i/£ is the pole nearest to the real axis and C is some constant. Thus we see 
that the freeze-out of the mode functions leads to the freeze-out of the correlation length 
£. The result of the numerical analysis is shown in figs, (l.g) and (2.b) which confirms 
this behaviour and provides the asymptotic value for £ « 1/h. From these figures it is also 
clear that the freeze-out time is given by the expansion time scale, 1/h. More precisely, the 



23 



numerical values for £ can be accurately reproduced by the following formula obtained by a 
numerical fit 

h£ ~ 1.02 + 0.2 In h + 0.06h + . . . . 



This situation must be contrasted with that in Minkowski space-time |2S| where the 
correlation length grows as £ (r) « yfr during the stage of spinodal growth. Eventually, this 
correlation length saturates to a fairly large value that is typically several times larger than 



the zero temperature correlation length ||28|| . We see that in the de Sitter case the domains 
are always horizon-sized. 

V. INFLATON ROLLING DOWN 

We now study the situation in which the inflaton zero mode rolls down the potential hill 
and consider initial conditions such that 7/(0) 7^ 0; r/(0) = 0. In this case there will be two 
competing effects. One will be the growth of spinodal fluctuations analyzed in the previous 
section, while the other will be the rolling of the zero mode. Which effect will dominate 
the dynamics is a matter of time scales. Before embarking on a numerical analysis of these 
cases, it is illuminating to try to understand under what conditions the dynamics will be 
driven by either the zero mode or the spinodal fluctuations. We will analyze first the case 
of the large N equations, and compare later to the Hartree case. 

The large N equations are summarized by equations Q3.12| ) - ( |3.14 ). The important part 



of the mode equations that determine the spinodal growth of long-wavelength modes is the 
term, if + gll. We can obtain an approximate estimate for the time scales on which the 
contribution of the zero mode becomes significant as follows. Let us consider the situation 
in which 77(0) << 1 and neglect the non-linear and backreaction terms (the last two terms) 
in the equation of motion for the zero mode (|3.12 ). At long times, but smaller than the time 



at which either the quantum fluctuation or 77 become of 0(1), we find, 

7](r) « r](0)e^- 3/2)hT } (5.1) 

with v given in eqn.( [4.3| ). The estimate for the time scale for the zero mode to be, rjf m 0(1), 
is approximately given by, 



(2i/-3)/i 



In 



Vf 
r/(0) 



(5.2) 



Comparing this time scale to the spinodal time scale given by ( J4.1U ) , for which quantum 
fluctuations grow to be of 0(1), we see that when, 
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V(0) « A 1/4 , 



(5.3) 



the quantum fluctuations will grow to be 0(1) much earlier than the zero mode for T > T c 
(for Tj = the bound becomes 77(0) << A 1 / 2 ). In this case the dynamics will be driven 
completely by the quantum fluctuations, as the zero mode will be rolling down the potential 
hill very slowly and will not grow enough to compete with the quantum fluctuations before 
the fluctuations grow to overcome the tree level terms in the equations of motion. In this 
case, as ar gued previously, the large N and Hartree approximations will be completely 
equivalent during the time scales of interest. 
On the other hand, if 

77(0) » A 1/4 , (5.4) 

then the zero mode will roll and become 0(1) before the fluctuations have enough time to 
grow to 0(1) (77(0) >> A 1 / 2 for Tj = ). In this case, the dynamics will be dominated 
by the rolling of the zero mode and is mostly classical. The quantum fluctuations remain 
perturbatively small throughout the inflationary stage which will end when the velocity of 
the zero mode modifies the equation of state to terminate de Sitter expansion. 

For r/(0) ~ A 1 / 4 (or 77(0) ~ A 1 / 2 for Tj = 0), both the rolling of the zero mode and the 
quantum fluctuations will give contributions of the same order to the dynamics. In this 
case, the quantum fluctuations will be large for the long-wavelength modes and the classical 
approximation to the inflationary dynamics will not be accurate. 

Since the scenario in which 77(0) >> A 1 / 4 , in which the dynamics is basically driven 
by the classical evolution of the zero mode has received a great deal of attention in the 
literature, we will not focus on this case, but instead analyze numerically the cases in which 
77(0) ^ but such that 7/(0) < A 1/4 . 



A. Numerical Analysis: 



We have evolved the set of equations of motion given by (|3.12|) and ( 3.13| ) numerically 



with initial conditions ( |3.14| ) for the large N case, and ( |3.15|) and ( |3.16| ), with the corre- 



sponding initial conditions ( |3.17| ) on the mode functions for the Hartree case. The numerical 



code is the same as in the previous section with the same relative errors. 
Large N: 

Fig. (3.a,b) show i](r) and g£(r) vs. r for the values A = 1CT 12 ; T/T c = 2; 77(0) = 
1CT 5 ; r)(0) = 0. Clearly the dynamics is dominated by the fluctuations; the zero mode 
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grows but is always negligibly small compared to gH(r). The time scale at which (?£(t) 
grows to be of order one is about the same as in the case, rj(0) = 0, and all the behavior for 
the mode functions, correlation length, energy density, pressure, etc. is similar to the case 
analyzed in the previous section. 

Asymptotically, we find that the sum rule, 

- 1 + r? 2 (oo) + #£(oo) = 0, (5.5) 

is satisfied to our numerical accuracy. This is the same as the situation in Minkowski space- 



time [ 25|)p6 |, and when rj ^ 0, this sum rule is nothing but the Ward identity associated 
with Goldstone's theorem. The fluctuations are Goldstone bosons, minimally coupled, and 
the symmetry is spontaneously broken with a very small expectation value for the order 
parameter as can be read off from figure (3. a). For r > r s , the dynamics freezes completely 
and the zero mode and the fluctuations achieve their asymptotic values much in the same 
way as in the case rj = studied in the previous section. Again, the correlation length 
becomes independent of time with £(oo) « 1/h in a time scale given by 1/h. 

Because there is a damping term in the zero mode equation, it is reasonable to assume 
that asymptotically there will be a solution with a constant value of n. Then the Ward- 
identity, 77(00) [—1 + r] 2 (oo) + <?£(oo)] = 0, must be fulfilled. In the large N case, the only 
stationary solutions are i) 77 = 0; g£(oo) = 1, or, ii) 77(00) 7^ 0; — 1 + ?7 2 (oo) + g£(oo) = 0. 
To have a consistent solution of the mode functions, it must be that the effective mass term 
[— 1 + ?7 2 (oo) + g£(oo)] vanishes asymptotically, leading to the mode equations for massless, 
minimally coupled modes which are asymptotically independent of time as shown in the 
previous section (see eq. (|4.24Q ). Furthermore, from fig. (3.b) it is clear that S'S(r) remains 



constant at long times, again unlike the case of free massless fields with Bunch-Davies 



boundary conditions in which case the fluctuation grows linearly in time (T^-|T4|. 
Hartree : 

Figure (3. a) also shows the evolution of the zero mode in the large N and Hartree case. 
Although there is a quantitative difference in the amplitude of the zero mode, in both cases 
it is extremely small and gives a negligible contribution to the dynamics. In the Hartree 
case, however, there is no equivalent of the large N sum rule; the only stationary solution 
for rj 7^ is, n 2 (oo) = l;<?£(oo) = 0. Such a solution leads to mode equations with a 
positive mass term and mode functions that vanish exponentially fast for r — ► 00 for all 
momenta. However, whether the asymptotic behavior of the Hartree solution is achieved 
within the interesting time scales is a matter of initial conditions. For example in fig. (3. a) 
the initial condition is such that the time scale for growth of the quantum fluctuations is 
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much shorter than the time scale for which the amplitude of the zero mode grows large 
and the non-linearities become important. In the large N case the sum rule is satisfied 
with a large value of the quantum fluctuations. In the Hartree case the equivalent sum rule 
— 1 + 3?yff + 3gY,H — is satisfied for a very small tjh and a ~ 1/3. The modes become 
effectively massless and they stop growing. 

The equation for the zero mode (see eqn. ( |3.15| )) still has an uncancelled piece of the 



non-linearity, — 2r/ 2 ; however the derivatives and the amplitude of rj are all extremely small 
and though the zero mode still evolves in time, it does so extremely slowly. In fact the 
Hartree curve in fig. (3. a) has an extremely small positive slope asymptotically, and while 
t)h grows very slowly, diminishes at the same rate. In the case shown in fig. (3. a), we 
find numerically that T]h/Vh ~ 10~ 7 at r = 150. Before this time most of the interesting 
dynamics that can be captured with a fixed de Sitter background had already taken place, 
and the backreaction of the fluctuations on the metric becomes substantial requiring an 
analysis that treats the scale factor dynamically. 

The conclusion of our analysis is that in the region of initial conditions for which the 
quantum fluctuations dominate the dynamics, that is for 77(0) << A 1 / 4 , both large N and 
Hartree give the same answer on the relevant time scales. The figures for H(r)/7i(0) are 
numerically indistinguishable from the case of figures (1). 

For comparison, we show in figure (4) both cases for the zero mode, for the same values 
of the parameters as in figures (1,3) but with the initial condition 77(0) = 10 -3 ; r)(0) = 0. 
This is a borderline case in which the time scales for the evolution of the zero mode and 
quantum fluctuations are of the same order and there is no clear separation of time scales 
between these two competing terms in the evolution equations. 

We see that in the large N case the zero mode rolls to a final amplitude which is 0(1) 
and of the same order as g£(oo) and the sum rule is satisfied. However, the Hartree case 
clearly shows the asymptotics analyzed above with ?7#(oo) = 1; gE#(oo) = 0. 

This particular borderline case is certainly not generic and would imply some fine tuning 
of initial conditions. Finally the case in which 77(0) » A 1 / 4 (or A 1 / 2 for T< = 0) is basically 
classical in that the dynamics is completely given by the classical rolling of the zero mode 
and the fluctuations are always perturbatively small. 

VI. DISCUSSION AND CONCLUSIONS 

We have identified analytically and numerically two distinct regimes for the dynamics 
determined by the initial condition on the expectation value of the zero mode of the inflaton. 
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1. When 77(0) << A 1 / 4 (or A 1 / 2 for Tj = 0), the dynamics is driven by quantum (and 
thermal) fluctuations. Spinodal instabilities grow and eventually compete with tree 
level terms at a time scale, t s > — 3/iln[A]/2. The growth of spinodal fluctuations 
translates into the growth of spatially correlated domains which attain a maximum 
correlation length (domain size) of the order of the horizon. For very weak coupling 
and h > 1 this time scale can easily accomodate enough e-folds for inflation to solve 
the flatness and horizon problems. The quantum fluctuations modify the equation of 
state dramatically and at this time scale can modify the dynamics of the scale factor 
and provide a means for a graceful exit to the inflationary stage without slow-roll. 

This non-perturbative description of the non-equilibrium effects in this regime in which 
quantum (and thermal) fluctuations are most important is borne out by both the large 
N and Hartree approximations. Thus our analysis provides a reliable understanding 
of the relevant non-perturbative, non-equilibrium effects of the fluctuations that have 
not been revealed before in this setting. 

These initial conditions are rather natural if the de Sitter era arises during a phase 
transition from a radiation dominated high temperature phase in local thermodynamic 
equilibrium, in which the order parameter and its time derivative vanish. 

2. When 77(C)) >> A 1 / 4 (or A 1 / 2 for T { = 0), the dynamics is driven solely by the classical 
evolution of the inflaton zero mode. The quantum and thermal fluctuations are always 
perturbatively small (after renormalization), and their contribution to the dynamics is 
negligible for weak couplings. The de Sitter era will end when the kinetic contribution 
to the energy becomes of the same order as the 'vacuum' term. This is the realm 
of the slow-roll analysis whose characteristics and consequences have been analyzed 
in the literature at length. These initial conditions, however, necessarily imply some 
initial state either with a biasing field that favors a non-zero initial expectation value, 
or that in the radiation dominated stage, prior to the phase transition, the state was 
strongly out of equilibrium with an expectation value of the zero mode different from 
zero. Although such a state cannot be ruled out and would naturally arise in chaotic 
scenarios, the description of the phase transition in this case requires further input on 
the nature of the state prior to the phase transition. 

We have learned from this work that non-equilibrium effects can alter scalar field dy- 
namics in a dramatic way, under well specified and physically reasonable conditions. In 
particular, this results bring up the tantalizing possibility that when the scale factor is 
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coupled to the inflaton dynamically, in a full backreaction treatment, some of the stan- 
dard results concerning the time evolution of the inflaton could be modified in unexpected 
ways. The fact that we found massless fields whose fluctuations did not grow linearly is 
an extremely interesting result, especially in light of how the quantum fluctuations become 
density perturbations. We are currently working on the formalism that allows us to couple 
gravity dynamically in a non-equilibrium way. Furthermore, it would be interesting to ex- 
tend the program of reconstruction of the inflaton effective potential (see |35[| ) to include 
the possibility of the dynamics being driven by spinodal fluctuations and not by "slow roll" 
of the inflaton zero mode. 

The results presented here raise some further interesting questions: how do these fluctu- 
ations contribute to the spectrum of primordial scalar density perturbations? How should 
the approach to cosmological perturbations based on slow-roll be modified to the case stud- 
ied here in which spinodal fluctuations, not slow-roll drive the dynamics? We are currently 
studying these issues within the consistent non-perturbative, non-equilibrium program pre- 
sented in this article. 
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Figure Captions: 
Fig. la: gE(r) vs. r for 77(0) = 77(0) =0; A = 10~ 12 ; r = 2; h = 2. 

Fig. lb: A [p(r) + e(r)] /2m^ vs. r for the same values of parameters as in fig. la. 

Fig. lc: H(t)/H(0) vs r for the same parameters as in fig. la. 

Fig. Id: H(t)/H 2 (t) vs r for the same values as in fig. la. 

Fig. le: gS(q,r) vs. q for r = 60. Same parameters as in fig. la. 

Fig. If: S(p,t)/S(0,t) vs. p for r > 2. Same parameters as in fig. la. 

Fig. lg: £(t) vs t. Same parameters as in fig. la. 

Fig. lh: In [|/ 9 (t)| 2 ] vs. r for q = 0,4, 10. Same parameters as in fig. la. 

Fig. 2a: #£(r) vs. r for 77(0) = 77(0) = 0; ; A = 10~ 12 r = 2; h = 0.1. 

Fig. 2b: £(t) vs r. Same parameters as in fig. 2a. 

Fig. 3a: rj(r) vs. r for A = 10~ 12 ; r = 2; /i = 2; 77(0) = 10~ 5 ; r)(0) = for large 

N (solid curve) and Hartree (dashed curve). 

Fig. 3b: (?S(r) vs. r for the same values of the parameters as in fig. 3a, for large N 
(solid curve) and Hartree (dashed curve). 

Fig. 4: Comparison of zero mode dynamics for large iV (solid curve) and Hartree 
(dashed curve) for 77(0) = 10~ 3 ; rj(0) = and all other parameters as in figures (1,2). 
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